Scale Free Cluster Distributions from Conserving Merging-Pragmentation Processes. 
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We propose a dynamical scheme for the combined processes of fragmentation and merging as a 
model system for cluster dynamics in nature and society displaying scale invariant properties. The 
clusters merge and fragment with rates proportional to their sizes, conserving the total mass. The 
total number of clusters grows continuously but the full time-dependent distribution can be rescaled 
over at least 15 decades onto a universal curve which we derive analytically. This curve includes a 
scale free solution with a scaling exponent of —3/2 for the cluster sizes. 
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The combination of the two basic dynamical processes, 
fragmentation and merging, is of importance in a number 
of situations of nature and society. Examples range from 
the dynamical evolution of companies where mergers and 
split-offs are important ingredients in daily business life. 
Empirical observations of US companies indicates that 
their sizes follow a scale invariant distribution An- 
other example comes from the dynamics of grain bound- 
aries in crystal growth where neighboring grains might 
merge at the same time as other grains fragment away 
from larger crystals. In crystal growth diffusion of grain 
boundaries can also be of importance in an interplay with 
merging and fragmentation |^ H, Fish schools are 
known to merge together resulting in scale free distribu- 
tion up to a school size after which it becomes exponen- 
tially distributed 0. In polymer physics merging and 
fragmentation are of importance in relation to gel and 
shattering transitions ; other examples are aerosols Q 
and scaler transport [g fragmentation/coagulation (i.e. 
merging) models (see e.g. j^, E| for reviews on merg- 
ing/fragmentation processes). Finally, we mention an 
astrophysical application for the formation of stars and 
interstellar clouds . 

It is known empirically that the resulting cluster dis- 
tribution in some of the above mentioned examples has 
scale free behavior. There has however, to the best of 
our knowledge, neither been proposed an analytical the- 
ory nor any fundamental model studies of these scale free 
distributions in closed systems driven by combined merg- 
ing/fragmentation processes. In this Letter we respond 
to this challenge by introducing a dynamical model sys- 
tem where clusters merge and fragment with rates pro- 
portional to the sizes, conserving the total mass of the 
system. Indeed, some previous models have been intro- 
duced to explain the scale free behavior but they relay on 
a source term for small cluster sizes and thus represent 
systems with non-conserved mass 0, 0, 0, 0, 0, . 
The dynamics of our model self-organizes the system to 
possess a stationary "continuous" scale free cluster size 
distribution over many decades and up to a characteristic 



size above which only the largest clusters will be situated. 
We introduce an analytical theory for the model based 
on Smoluchowski's continuous cluster fragmentation and 
coagulation (i.e. merging ) model 113. We show that this 
model with "mass" kernels exhibits an exact scale invari- 
ant solution of the cluster distribution, n{x) ~ x~^^^, 
at a critical point where merging processes balance the 
fragmentation processes. In this solution, it turns out 
that there is a gap in the distribution above which only 
a few large clusters exist. The scaling exponent — | is in 
agreement with distributions of fish schools and many 
critical branching processes ^3|> foi' instance. 

Our model is defined in terms of M clusters 
each of size Xi, where tilde ("') in the following is used 
for dimensionfull quantities. The total initial mass of all 
clusters is thus M = Y^'^i ^i- time i the clusters i 
and j ^ i are chosen to merge with a rate Kij — (3xiXj 
and the cluster k is chosen to fragment with the rate 
Fk — fxk- The size of the two clusters resulting from 
a fragmentation is chosen randomly (ie. with a uniform 
distribution). The average time 6t for the next merging 
or fragmentation event to occur is obtained from the total 
merging and fragmentation rates, fcmerge and k^j-g^g, as 

Si = (%ag + fcmerge)~\ where 



fcmerge ^ 1^1^ (^^' " 



and 



%ag = E^' = /^ 



and M2 — is second moment of the distri- 

bution. In time-step 5t the merging and fragmenta- 
tion processes are thus associated with the probabilities 
Pij = KijSt and pk = FkSt, respectively. Since the pro- 
cesses preserve the total mass we can express the dynam- 
ics in terms of the dimensionless quantities, 
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where b = M/S/f and Ah = J2i ^l- 
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In Fig. ^ we show the cumulative distribution, 
N{x,t) = Q{xi{t) — x), for a simulation with 6 = 3 
at various times t. Here, Q{x) is the Heaviside func- 
tion, Q{x) = 1 for X > and Q{x) — for x < 0. 
At each time t, the fragments larger than a characteris- 
tic cross-over scale, x ^ Xc(t), are distributed according 
to a power law, N{x) ~ x" with a = —1/2. As seen 
in the figure, xdt) for t ^ oo, implying the exis- 
tence of a scale invariant state which extends to smaller 
and smaller scales as time increases. The same data are 
shown in Fig. |21 with N{x,t) as function of t for vari- 
ous values of x. Interestingly, the total number of do- 
mains M{t) — iV(0, t) grows as A/'(t) ~ with the power 
7 = —a ~ 1/2. For each a; > 0, the distribution is 
time-dependent, N{x,t) ~ t^^"^, for times smaller than 
a characteristic cross-over time t <C tc{x), and becomes 
time independent, N{x, t) — N{x) for t 3> tc{x). 

Simulations for other values of 6 > 2 lead to distri- 
butions, N{x), and evolutions Af{t) which are very close 
to power-laws, N{x) - and 7V(t) V'-^'^ with b- 

dependent exponents 7(6) and a{b) satisfying 7(6) = 
—a{b) 0. Nonetheless, only at 6 = 3 does the scale- 
invariant behavior become exact. For 6 < 2, the maxi- 
mum mass mi approaches mi = as i — > 00 |19|. im- 
plying that no time-independent solution exists in this 
parameter region. 

To analyze the origin of the scale invariant solution we 
express the dynamics in terms of Smoluchowski's contin- 
uous cluster fragmentation and coagulation (i.e. merg- 
ing ) model [13. We have previously presented this type 
of rate equations including diffusion, fragmentation and 
coagulation (ie. merging) terms, which are of relevance 
for dynamical processes of ice crystals, alpha helices and 
macromolecules mill. Let n{x,t) = -d^N{x,t) de- 
note the number density of clusters of size x at time 
t. Smoluchowski's fragmentation-coagulation model de- 
scribes the dynamics of clusters subject to fragmentation 
and coagulation processes. If K(x,x') and F{x,x') rep- 
resent, respectively, the rate of coagulating two clusters 
of size X and x' into a cluster of size x -\- x' and the rate 
of fragmenting a cluster of size x + x' into two clusters of 
size X and x' the model takes the form 



dtn{x,t) = —n{x,t) j F{x — x' ,x')dx' 



2 / F{x,x')n{x + x',t)dx' 

1 r 

— \ K{x — x' , x')n{x — x' ,t)n{x' ,t)dx' 

2 Jo 

n{x,t) [ K{x,x')n{x',t)dx', (2) 



where L represents the largest possible cluster size. Using 
dimensionless quantities, Eq. our model corresponds 
to the choice, L = 1, F{x, x') — 1 and K{x, x') — b-{x-x'). 
The general form of the merging-fragmentation process 
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FIG. 1: Cumulative distribution, N{x,t), of the fragment 
sizes X at various times, t = 5-10', 2-10*, lO'', 5-10^, 2.5-10^°. 
The thick solid line is the average cumulative distribution, 
N{x), for X > 10~® taken over times t > 10^. At any instant, 
A^(a;, t) displays a characteristic cross-over from an almost flat 
distribution to a scale-invariant distribution, N{x) ~ x~^^'^, 
at larger x. As time goes the scale free part of the distribu- 
tion extends to smaller and smaller scales. The inset shows 
the behavior of N{x) in the large part of the mass spectrum. 
The dotted line is the simple scaling form, Ni(x) = Ax~^^^, 
with A ~ 0.9, and the dashed line is the predicted cumu- 
lative distribution, No{x) = A (x^^''^ — S/lBAvr) (see text). 
The dots in the figure represent the typical size (50% frac- 
tile) of the 10 largest fragments, ie. points corresponding to 
N{mi) = i — 1/2, where rrii the the typical size of the i'th 
largest fragment. The predicted cumulative distribution, A'o, 
deviates from data around x = ? ~ 0.1 corresponding to the 
typical size of the third largest fragment, ms « 0.09. The 
naive scaling form, A^i, deviates from data already around 
X ~ 0.02. 
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FIG. 2: Cumulative distribution, N{x,t), as function of time 
t at various sizes, x = 0, 10"^°, 10"^, 10"*, 10"^, 10"^. For 
each positive value of x, N{x,i) displays a cross-over from 
a power law N{x,t) ~ t^''^ at short times to a constant, 
N{x,t) = N{x), at large times. 



3 



then reads 



dtn{x,t) = —xn{x,t) + 2N{x,t) 



+ - / x' {x — x')n{x' ,t)n{x — x' ,t)dx' 
2 Jo 

— bxn{x,t), (3) 

where N{x,t) = 72{x' ,t)dx' and the mass normahza- 

tion, xn{x,t)dx — 1, has been used. We note that 
the representation of the dynamics in terms of a density 
function, n{x,t), is only accurate if the typical spacing Sx 
between clusters of size ~ x, satisfies Sx <^ x, implying 
X ^ 1. In particular, the assumption of smoothness is 
important for the integral expressions for the fragment 
gain and loss due to merging to be correct. These terms 
erroneously include 'self-merging' processes, x + x ^ 2x, 
as well, thereby giving an addition to the actual merg- 
ing rates which can only be ignored in the limit 5x <C x. 
Eq. (|2I3|I are therefore not adequate for the dynamics of 
fragments carrying a large fraction of the total mass, as 
the error in the merging rates become non-negligible here 
and fail to guarantee total mass conservation. In fact, 
Eq. ©J is only mass conserving in the limit L — > oo, for 
which the particular rescaling to dimensionless quanti- 
ties in Eq. J^l, can not be done. In the following, we will 
tentatively assume the continuous representation to be 
correct up to some length scale I <C 1 and seek analytical 
solutions to Eq. 10) in this regime. 

Critical points are generally associated with scale in- 
variant behavior. For cluster dynamics with merging and 
fragmentation, scale invariance has been obtained previ- 
ously in [H 111 El lisl ll^ . although all under some- 
what different circumstances i.e. with source terms and 
thus not in a closed system. Therefore, we are search- 
ing for stationary scale-invariant solutions on the form, 
no(a;) ~ a;"^^ or Nq{x) = Ax°' + C. Inserting this scaling 
relation into Eq. ||2J) we obtain 

= -Ax" {a + 2 + ba) + 2C 

+ ^(aA)2x(2"+i)*S(a + l,a-Hl) (4) 

where B{a+1, a+1) is the Beta function B{a+1, a+1) — 
T{a + l)T{a + l)/r(2Q; -I- 2) under the assumption that 
a+1 > 0. Power counting arguments now easily gives 
that only two values for the exponent is possible, namely 
a = —I or a = —1/2. The value a = —1, however will 
ruin the assumption a -I- 1 > and is inconsistent with 
the finite mass of the system |23|. Thus we come to the 
conclusion that there is only one scale invariant solution 
to Eq. namely no{x) = A/2x~^^^, which becomes 
an exact solution provided that a + 2 + ha — or b — 3. 
Also, stationarity requires C = —3/16A^Tr so the total 
number of clusters larger than x for Xc{t) <^ x < I, is 
simply 



The actual value of A depends on how much of the to- 
tal mass is situated in the discrete part of the spectrum 

given by the few fragments with x > I, J^^ xnQ{x)dx + 
>i 2^4 — 1- Setting (incorrectly) I — 1 implies A — 1. 
To assess the scaling behavior in the large mass limit, we 
show in the inset of Fig. ^ the average cumulative dis- 
tribution for X > 10~^ where size fluctuations becomes 
noticeable. Using a fitted value of A w 0.9, the analyt- 
ical solution, Eq. ||SJ) is observed to be accurate all the 
way up to Z ^ 0.1, as expected from the condition I <^ 1. 
The value for / corresponds to the typical size (50% trac- 
tile) of the third largest fragment, « 0.09, effectively 
meaning that only the two largest fragments belong to 
the discrete spectrum. In fact, the distribution N{x) has 
significant gaps around the two largest fragments, e.g. 
6x = l/2(mi - mg) « 1/2(0.47 - 0.09) = 0.19, com- 
pared to the small differences Sx ^ 1 around all other 
fragments. 

Eq. lO suggests that the full time dependent distri- 
bution, N(x,t) can be expressed in the form N{x,t) ~ 
(j){xt) + C. Indeed, setting z = xt one obtains from Eq. 

@ 

-tmj + i)cf>\z)+zr{z))^ 

{z(t,'{z) + 20(z) + bz(l)'{z)) + 
C+l/2t^->-^b j\'{z- z')c^'{z)c^'{z- z')dz\ (6) 

JO 

which is a function of z only, provided that either 7 = 1 
or 7 = 1/2 and l/2b£z'{z - z')(l,' {z)(j)' {z - z')dz' « 
—C. Requiring N{x,t) — No{x) for x ^ Xc{t), implies 
(j){z) ~ z" for large z and 7 = —a — 1/2. Neglecting the 
two last terms on the rhs. of Eq. one can obtain 
the distribution for all z by solving the linear differential 
equation, 

-3/2(/)'(z) - z(j)"{z) = z^'{z) + 2(/)(z) + 3z(j)'{z), 

where 7=1/2 and & = 3 has been inserted. The solution 
yields (/)(z) = aiz^^^'^ + a2ip{z), where V'(^) = 
and erf(y) is the error function, erf(?/) = e^*" dr. 

The integration constants, ai and 02, can be found from 
the known behavior of 4>{z) in the limits z ^ 00 and z — > 
0. Requiring 0(z) — s- const for z — > implies ai = and 
since 'ip{z) — > z~^/^ for large z, one obtains 02 = A. For 
this solution the convolution term in Eq. is bounded, 
< !^ z'{z - z')(j)'{z)(j)'{z - z')dz < 5(1/2, 1/2) and has 
a weak z-dependence which is negligible at larger times. 
The final form of N{x,t) then reads 



N{x, t) = t^l'^4>{xt) - 3/16^V, (t){z) = A- 



rf(2^/i) 



No{x) = A {x-^/"^ - 3/16A7rj . 



(5) 



(7) 

In FigureEl we show {N{x,t) + 3/16A^Tr)t-^/^ , as func- 
tion of z — xt for 50 distributions taken in the time 
interval 10^ < t < 2.5 • 10^°. The thick solid fine is the 
predicted functional form, 0(z) which agrees perfectly 
with data. The full time dependent solution, Eq. {Tj), of 
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FIG. 3: The rescaled cumulative distribution, t~^^'^(N{x,t)~ 
C), where C — —3/16A^-k, Eq. Q, as function oi z = xt 
taken at 50 different times in the interval 10^ < t < 2.5 ■ 
10^°. In black solid line is shown the predicted functional 
form0(^)=vl2rf{^. 



the combined process of fragmentation and coagulation is 
the main result of this Letter. The spread around (f>{xt) 
for large values of xt is caused by noticeable fluctuations 
of the largest fragments of the distribution. Eq. 
also shows that the characteristic cross-over scale, Xc{t) 
(above which the the system is stationary and scale free), 
goes with time as Xc{t) ^ t^^. Conversely, had we intro- 
duced a smallest length scale, 6, in the system, we would 
expect convergence to a full stationary distribution after 



time t ^ S ^, displaying scale invariance in the region 
(5 < X < 1. 

Let us now comment on conserved versus non- 
conserved merging/fragmentation processes. Imagine a 
simulation where clusters, when becoming too large, 
above the scale I, are disregarded and the total mass is 
not conserved. From these kind of simulations we have 
observed a transient scale invariant cluster size distribu- 
tion, n{x) ~ x~'^ , consistent with simple power count- 
ing of Eq. |(2Jl (corresponding to a = —1). However, 
the total mass will eventually become smaller than I, af- 
ter which the dynamics again is mass conserving leading 
to the n{x) ^ x~'^^'^ distribution. The —2 versus —3/2 
solutions might be relevant in discussions of data from 
a variety of systems depending on the degree to which 
they are mass conserving. In the case of American com- 
panies, the exponent of the scale invariant distribution 
is closer to —2 than —3/2 [J, which in the light of our 
model indicates that the total wealth of companies is not 
conserved. 

In this letter we have presented a general formalism 
for scale invariant and critical behavior in systems where 
merging balances fragmentation. We have proposed a 
dynamical model where the merging/fragmentation pro- 
cesses drive the system into a critical state. From an ana- 
lytical solution of Smoluchowski's equation we show that 
this scale invariant state represents a unique part of the 
cluster distribution. We suggest that the proposed prin- 
ciple could be the dynamical background behind many 
scale free systems of society and nature. 

We are grateful to S. Manrubia, S. Maslov and K. 
Sneppen for useful discussions. 
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To be precise, a — ~1 implies that the total mass 
in the scale-invariant region, Xc{t) < x < I, becomes 
f^^^xno{x)dx ~ \og{l/xc{t)) — > CX3 as the lower scale, 
Xc(t), tends to zero at large times. 



